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DESCRIPTION OF COURSE 

MPE 251: Metallurgical Thermodynamics : 3 credits 

Principles of Thermodynamics; The laws of thermodynamics; thermochemistry, enthalpy and 
heat balances, standard state; 

Chemical equilibria, entropy and equilibrium state, equilibrium constant, free energy change 
in chemical reactions; Chemical equilibrium and metal extraction from aqueous media (Eh- 
pH or Pourbaix diagrams); Chemical equilibrium and stability of metal oxides (Ellingham 
diagrams) 

Solution Thermodynamics ; gas mixtures; aqueous solutions; the structure of solutions; dilute 
and concentrated solutions; the activity; the Gibbs-Duhem equation; Debye -Hiickle theory, 
solubility; complex formation; solid solutions and melts; 

Phase Equilibria, fundamentals of phase diagrams; 

Electrochemistry, driving force of electrochemical reactions, electrode potentials and 
electrochemical series, the effect of concentration on the EMF, galvanic and concentration 
cells, polarization and overpotential; electron activity and introduction to p£-pH diagrams. 
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1 INTRODUCTION (BASIC CONCEPTS) 


1.1 System 

A system is any object, quantity of matter, region, etc., selected for study and mentally set apart from 
every thing else (surroundings). The system and its surroundings are separated by the boundary. The 
boundary serves to isolate the system from its surroundings or to provide for interaction in specific 
ways between the system and the surroundings. 

An isolated system has no interactions of any kind, neither exchange of matter nor energy with its 
surroundings. When there is exchange of matter and energy the system is said to be an open system. 
and when there is exchange of energy but no matter, then the system is a closed system . 

A homogeneous system is one in which every property of the system is constant irrespective of the 
coordinates of the system. In a non- homogeneous system at least one property has a continuous 
function with the coordinates (place). A heterogeneous system is a system in which at least one 
property is a discontinuous function with the coordinates. 




p la ce 

Horn o g e n e o u s 




p la ce 

Heterogeneou 


1.2 The Equation of State 

When constraints are laid upon a system, one can find a unique relationship between the 
thermodynamic variables, eg. for a given mass of oxygen there is a unique relationship between the 
pressure, volume and any other thermodynamic variable. This relationship is expressed by an 
empirical equation called the Equation for Functionl of State. 

f(p r V,x)= 0 x is a third variable e.g. temperature 

The function of state describes the condition of a system as it is now, without reference to prehistory. 
Therefore the change of state (from state 1 state 2) is independent of the path, but only on the 
initial and final stages. 

A simple example of an equation of state is the ideal gas law: pV — nRT 

A system has an identifiable, reproducible state when its properties are fixed, or in other words the 
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properties of a system are fixed by its state. 

1.3 Extensive and Intensive Properties 

Certain properties (or variables) of a system show dependence on the size of the system and are 
proportional to the mass of the system. These properties are termed extensive properties . These 
properties are additive eg., number of moles of the system (N), the internal energy (U), the volume 
(V), the entropy (S), etc. 

On the other hand variables such as temperature (T), pressure (P) are independent of the size or mass 
of the system. These properties are called intensive properties. 

1.4 Equilibrium 

Any change in any system can be measured by determining changes in the property of the system, e.g 
by using such instruments as thermometers, pressure gauges, etc. The concept of equilibrium is 
central in thermodynamics. When a system is displaced from a state of equilibrium, the system 
undergoes a process or series of processes during which its properties change until a new equilibrium 
is reached. A system is in equilibrium when none of its thermodynamic properties changes with time 
at a measurable rate. During such a process the system may interact with its surroundings so as to 
interchange energy in the form of heat and work . Such changes do not continue indefinitely, but come 
to a state of internal equilibrium . If the changes are due to external interactions, then external , as well 
as internal equilibrium is reached after some time. The equilibrium state depends on the constraints 
(conditions) imposed on it and generally there is only one true equilibrium state for a set of 
constraints. 

1.5 The Concept of Work 

In Physics the definition of work is: 

work done (*Wi) = force x distance 
ie: 5wi=Fidxi = PidAidx, = -pidV, 


Total work done is: p 1 = p p' > p 

5w = -p'dV 

w 

For a measurable volume, the total work done on a system is given by: 

Note: *w indicates an inexact differential (that is it is not the differential of any function). In other 
words, it means that the integral depends on the path taken, eg. heat associated with work. 

In general, if an expression, df = M(x,y)dx + N(x,y)dy represents a differential of some function f, then 
it is called an exact differential . An inexact differential can be converted into an exact differential by 
multiplying by some integrating factor e.g. 2ydx + xdy is inexact because there is no function which 
satisfies the condition: 

d M ( d N , However, if this expression is multiplied by x the resultant expression becomes 
v g / — v q x I an exact differential. 




2xyd x+x 2 dy 


becomes an exact differential of the function "x 2 y 
















Examples of work: 

1) Electrical work done in charging a capacitor is: 


Q2 

w= J i? dq 

Qi 


where q is the charge and 


' is the applied potential difference. 


2) Work done on a point of charge (q) moving through an electrical field (E) is: 

bw=F.ds=qE.ds =—qd S* 

3) Work done in a gravitational field is: 

2 

j {mg) 

p 'dV = - - Adh-mg 


1.6 Reversible and Irreversible Processes 

A process is reversible if the system remains infinitesimally close to equilibrium at all times with 
changes taking place at near zero rates. The system thus passes through continuous sequence of 
equilibrium states and can be reversed in the opposite direction along the same path without any input 
of additional work above that for the forward process. At any point in the process any intensive 
property (p,T) is continuous across the boundary. An example is the expansion of gas through volume 
Vi to V 2 assuming frictionless piston and that the applied pressure is the same as the pressure in the 
system. Reversible processes produce maximum work. Another example is tug-of-war with equal 
opponents. 

In irreversible processes the system passes through non-equilibrium states with unbalanced forces 
giving the process a unique direction (eg. the flow of water to a lower level, or the flow of heat to a 
colder body, or the dissolution of salts, or the reaction between ammonia and hydrochloric acid). For 
an irreversible gas expansion the applied pressure is not the same as the pressure in the system. 

Real processes have measurable rates of change, but because most real processes have their 
characteristics (eg. values of work, heat evolved) close to those of reversible processes, we assume 
reversible conditions for these processes for the purposes of calculation. 


2 

J dh = mg A h 



2 THE LAWS OF THERMODYNAMICS 
2.1 Zeroth Law 

The Concept of Temperature : The basis for empirical temperature scales is given by the so-called 
Zeroth Law of Thermodynamics, which states that: "Two bodies, each in thermal equilibrium with a 
third body, are in thermal equilibrium with each other. These bodies have a property in common, 
defined as temperature (degree of hotness." Different temperature scales have been developed (eg. 
Fahrenheit, Centigrade or Celsius) but even with these scales the measured temperature depends on 
the thermometer used (Hg-in-glass, or thermocouple, etc.). The Kelvin scale based on the ideal gas 
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law, the Second Law of Thermodynamics and the Kinetic Theory, is the most accurate one and is the 
scale used in thermodynamic calculations. 

2.2 Postulates 

The Laws of Thermodynamics constitute varieties of observations, abstract principles which are 
ultimately justified by the agreement of conclusions drawn from them with experiments. To explain 
these observations scientifically a number of postulates have been advanced: 

1) that there exists a form of energy, internal energy (U) which for systems in equilibrium state is an 
intrinsic, ("belonging naturally") property of the system, functionally related to the measurable 
coordinates which characterize the system, ie. 


U = U (S,V,Nj) where j = 1, 2, 3,..., k (for an open system). 


U = U (S,V,) for a closed system 

In the case of work being done on an adiabatically contained body of constant potential and kinetic 
energy, the pertinent function which describes the state of the body is the internal energy (U). 

Internal energy is a function of state and therefore any change in the internal energy is independent of 
the path, but only on the initial and final states. The internal energy is different from the kinetic or 
potential energy, which are external forms of energy. 

2) that the total energy of any system and its surroundings is conserved ( First Law of 
Thermodynamics ). In other words, the expenditure of a given amount of work, no matter what 
origin, always produces the same quantity of heat. 

For an isolated system the first law requires that its energy be constant. For a closed system the 
first law requires that energy changes in the system be exactly compensated by energy changes in its 
surroundings in the form of work and heat. (This law serves as the basis for the preparation of 
material and energy balances). 

Mathematically, for a closed system of constant mass in which only the internal energy changes: 

dU =bQ—b W where *Q represents the heat absorbed by the system from its 

surroundings, and *W is work done by the system. After integration the 

expression becomes: 


A U=Q-W 


For an isolated system change in internal energy is a result of "adiabatic" work done in bringing the 
system from initial to final equilibrium state, ie: 

A u=U-U=W ad 

If we consider all the different kinds of work (volume, surface, electrical, etc.) that can be done by the 
system then the first law can be written as: 

dU — bQ — pdV+y dA — IT dq —... + ^ dNj 

3) There exists a property, entropy (S), which for systems in equilibrium is an intrinsic property and 
is functionally related to the measurable coordinates of the system. 



S=S[U, V,Nj) 

According to Gibbs, entropy is the measure of the mixed-up-ness. It is a function of state. For 
reversible processes entropy change is calculated along the reversible path as: 


dS= 


8 g 

T 


T is the absolute temperature, and Q rev is the heat evolved or 
absorbed by the system. 


4) The entropy change of a system and its surroundings, considered together, resulting from any real 
process is positive and approaches a limiting value of zero for any process that approaches 
reversibility (The Second Law of Thermodynamics 1 


( or "the complete conversion of heat into work is impossible without causing some effect to occur 
elsewhere" or " no engine, operating in cycles, can absorb heat from a reservoir at uniform 
temperature and convert it into an equivalent amount of work without causing some change elsewhere 
in the universe") 


The Second Law deals with the conversion of heat into work and the law indicates that such 
conversion cannot be carried out to the full. For an isolated system the Law requires that the entropy 
increases or in the limit remains constant when the system is in equilibrium. For a closed system the 
Law requires that any entropy decrease in the system or its surroundings be compensated in the other 
part or in the limit remain constant at equilibrium. 

Every substance has a finite positive entropy, but at the absolute of zero of temperature the entropy 
may become zero and does so in the case of a perfectly crystalline substance - The Third Law of 
Thermodynamics . 


5) The macroscopic properties of a homogeneous (PVT) system in equilibrium can be expressed as 
functions of P, T, N only. 

2.3 Mathematical formulations of thermodynamic functions 

The most common systems encountered in the application of the 1 st and 2 nd Laws are gasses and 
liquids. The variables most often used are p, V, T and N (composition) and the work done in reversible 
processes in such systems are: 


bW re =pdV 


For closed systems: dU=bQ iev —bW iev or dU = TdS — pdV 

For open systems: U = U(S, V, Nj) 


dU = 


dUp 

dS' 


dS+ 


V,N, 


dU 

dV 


dV+Ya 


dU 


S,N , 


dN 


dN, 


j S.V.N 


By definition, 


I™.) =T; 

d S' 

oo v,n 


,du_ 

'dV 


=~p: 


S,N , 


.dU 

'dN : 




S.V.N. 


where p is the chemical potential 


p = f (T,P, phase) 
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The chemical potential plays a vital role in phase and chemical equilibrium considerations and is 
generally defined as the free energy per mole of a pure substance. 


The change in internal energy becomes: 

dU - TdS-pdV + X fly dNj This equation is known as the Fundamental 

Thermodynamic Equation (Relation) which on 
rearranging gives: 

dS=-dU+ — dV-^ J -dN, 

T T T 


showing S = f(U, V, Nj), which on differentiaton gives: 


dS= 


H) rfr-I' 


dS 


d U V,N. d V U,N 


dS 

dN : 


Comparing the last two equations it is clear that: 


) dN, 

U ,V ,N. 


_1 idS\ _P / 35 \ _ _ h/ 

W } UiN ~T : [ dN/ VJiN ~~F 

In terms of mole fractions (x), the Fundamental Thermodynamic Equation becomes: 
dU — TdS—pdV + ^ Hjdxj 


Formulations with other thermodynamic functions include: 

1) Enthalpy (H): H = U + pV dH=dU+pdV+Vdp 

Substituting dU from the Fundamental Thermodynamic Equation ans simplifying gives : 
dH = TdS + Vdp + ^ u y dN j 

from which it is evident that enthalpy is a function of the entropy, pressure and number of moles of the 
system, or 

H = H (S, p, Nj) (students are to find the relationship between enthalpy and pressure, 

temperature and the number of moles of a system, using the example 
shown with respect to entropy above) 

2) The Helmholz function (A): A = U - TS dA =dU — TdS—SdT 

Using similar substitution from the Fundamental Thermodynamic Equation, 

dA+ — SdT — pdV + ^ u ; dN j 



A = A (V, T, Nj) 


3) The Gibbs function (G): G = U + pV - TS = H - TS 

With similar manipulations as above, 

dG = Vdp - SdT+Y ii ; dNj G = G(p, T, N) 

This equation shows that the Gibbs energy is a function of the pressure, temperature and the number 
of moles of the system. 

All these energy functions are important in the study of thermochemistry, among others. For example, 

(i) it is easier to calculate the enthalpy change of an isothermal, isochoric system using the 
expression for (dH) if the entropy, pressure and number of moles of the system are changing. 

(ii) for an adiabatic, isobaric system the Helmholz function becomes important. 

(iii) the Gibbs function provides the criterion for deciding whether or not any kind of change 
(chemical, phase, etc) is possible, 

AG < 0, indicates a possible spontaneous change 

AG = 0, shows an equilibrium situation 

AG > 0, indicates that change is not possible, but rather the reverse change 
may tend to occur. 

2.4 Maxwell Relations 

These show the relationship between different thermodynamic expressions. These make it possible to 
substitute one expression for the other and make thermodynamic calculations easier. 

For example, from the expression: 


H = U + pV, 


dH — TdS + Vdp +X \ij dN j 


Also from H = H (S, P, Nj) dH = (^Mr) dS+(^~) dp + Y 


dS 


Therefore 


(^) =T; (^) -V; 


dp s. Nj 

,dH 


dH \ mu 


dS 


p,N, 


dp 


S.N, 


N ’ rj 

Uly j S.p.N 


Using similar manipulations with the Helmholz function, 

A = A (T, V, Nj) and dA = - SdT - pdV +X P 7 dN j 


dA 

dT 



,dA 

W 


=-p; 


T.N. 


, dA 


=!T, 


T.V.N . 


From G = G (T, p, Nj) and 


dG = Vdp— SdT + X P, dNj 
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dG 

dT 



Th above shows that, 


,dG_ ] 

y dp 


= V; 


dG 


T,N 


y dN. 


=h, 


; t.p.n, 


,dU 

'dS 


V.N. 


, dH 

[ dP 


s,N, 


t d_H_ \ ,d(T 

dp s.Nj dp TN . 



d_G_\ 
d T p 


2.5 Applications of energy functions 

The energy functions are introduced as a matter of convenience. We are mostly interested in, for 
example, the direction of changes between phases, chemical reaction possibilities, phase equilibria, 
diffusion, etc., as well as to derive important relationships between properties of substances. These 
functions will enable us to express the equilibrium conditions for example of these changes, among 
others, and to work with the functions which are connected to the process under study. 


2.5.1 Expressing U, S as functions of T and V 

From postulate (5), at constant composition (N = constant) the molar thermodynamic properties are 
functions of T and P. However, sometimes it is more convenient to use T and V as the independent 
variables instead of T and P. 

Using U = U(T,V) and S = S(T,V), after differentiation: 


dU = 


dU 

dT 


dT + 


dU 

dV 


dV 


dS= 


d_S_> 

dT' 


dT +( 


dS 

dV 


) dV 


By definition the molar heat capacity at constant volume (C„) is the temperature derivative of the heat 
added to a mole of material at constant volume and composition and it is a property of the material 
and a function of T, P, and N: 


r —( ^ ^ rev \ 

v ~ [ dT’ v ~ [ dT V 

From dU = TdS - pdV, 


dU 

dT 


) =T 


dS ' 

dT’y 


and 


dU 

dV 


=r 


ds 
d V ' 


From Maxwell relations, however, 


Therefore, 


dS 

dV 



V 


dU 

dV 


) =T 



or 


dU — C v dT + 



) ~p]dV 


Similarly, 


d p 

ar 


. (dl) 

“wy ar t 

dS = ~~ dT +(~~) dV 

T y dT v 


and therefore 



2.5.2 H, S as functions of T, p 

Similarly at constant N, 

H = H(T, p) and S = S(T, p) 

Similar manipulations as above can be used to find the relationship between enthalpy (H) and 
temperature and pressure, and the relationship between entropy (S) and temperature and pressure 
(students may want to do that on their own): 

dH = C p dT+[V-T{]dp 


C p = f(T, p, Nj) and by definition 

So we see that changes in entropy (dS), enthalpy (dH) and internal energy (dU), which are not 
measurable quantities can be expressed as measurable quantities. The pressure or volume dependence 
of molar heat capacities ( Cp and Cv) are determined from PVT data. The temperature dependence is 
given empirically as: 

C p = a + bT + cT 2 + dT 3 +... 

C p = a + bT + cT' 2 +... 

As T = OK C p vanishes, rising to a plateau at high temperatures with discontinuities occurring at the 
melting and boiling points. 

In a homogeneous phase at constant volume : 

Heat added to the system is given by T l 

q — A U = J C v dT 

Tl 


dS=^dT -{|£) dp 

_/ 6 H \ _/ ^ Q rev | 

dT dT 1 


Entropy change is, 


T2 

A Sy=i 


Tl 


Cy 

T 


dT 


At constant pressure 

Heat added to the system is given by: 


T2 

q = AU=j C p dT 

Tl 


Entropy change is 


T2 

A S= \ — E dT 

P J r r 

Tl 1 


At constant temperature 
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pi 

A G r — J Vdp 

p2 

At constant entropy (adiabatic 


A S t =-R In — 
Pi 


AS t =R In 


process') 


V2 

V l 


V , 

A U S --RT In—=■; 

V, 


A H,=RT In — 


2.6 The concept of Standard States 

Standard states are reference points to which energy functions can be related. Since we often consider 
changes in state, it would be absurd to consider "absolute values of these quantities. Usually, standard 
values are given at 101.3 kPa (one atmosphere) pressure and at 25°C (298.15 K). Each element in its 
most stable form at 25°C and 101.3 kPa is assigned a quantity of zero. The standard values are 
tabulated in Thermodynamic Table of Constants for substances rather than reactions. 

Other applications of the First Law are found in calculations in Thermochemistry, ie. heats of reaction, 
combustion, neutralisation, formation, decomposition, sublimation, vaporization, bond energies, 
crystallization, dilution, dehydration, among others. 


Examples of calculations in Thermochemistry 
.Material and Energy Balance calculations 





3 CHEMICAL EQUILIBRIUM 

Chemical equilibrium can be considered as a consequence of the Second Law of Thermodynamics. For an 
isolated system spontaneous, irreversible processes are accompanied by entropy increase in the system. 
This increase will persist so far as spontaneous changes occur, until no further spontaneous changes are 
possible, when the entropy increase will have reached a maximum. Then no further increase in entropy will 
be possible. 

In this case dSu.v = 0. 

The internal energy change in such a system will be: 

dU = TdS - TdS' - *w (TdS' is a contribution from the 
spontaneous changes). 

At reversible conditions TdS' = 0, (ie. dU = TdS - *w = TdS - pdV) 
for a closed system. A criterion for equilibrium therefore is: 

dU(S,V) = 0. 

Similarly other criteria for equilibrium are: 
dH(S,P) = 0; dA(T,V) = 0 

Another approach to chemical equilibrium is through the Gibbs free energy change. The trend of a reaction 
or change will always be in the direction for which change in Gibbs free energy (AG) is negative, and when 
equilibrium is reached the AG is zero. 

Consider the general equation: 

aA + bB = cC + dD where a, b, c, and d are stoichiometric 
coefficients. A,B,C,D are the various species. 


From dG = Vdp - SdT + X P, dN j 

the change in free energy at constant T and p is: dG = 'Y. u , dN . 

In an infinitely large system where finite changes occur without changes in the chemical potentials, for each 
mole of A consumed, 

A G=c\i c +d p^-up^-^p^ 

But n j = n* j+ RT \n(a j) 

Therefore substituting and simplifying gives: A G =A G°+RT In Q 
Q is the activity quotient Q = 


a b 



At equilibrium (for a closed system) AG T>P = 0 , or 

A G °= — R Tin Q e = — Rt In K where K is the equilibrium constant 
The equilibrium constant provides mathematical tool for theoretical predictions and likelihood of reactions. 
In general, for the reaction: 


aA + bB + ... = pP + qQ + 


A G°——RT\n- 


Up . Ur 


a A • ci B 


from which 


A G° = — RT In K„ 


where K a is the equilibrium constant expressed in terms of 
the activities and is a dimensionless ratio. 


The above equation is termed Reaction Isotherm Equation or the van't Hoff ISOTHERM. Ka is the 
equilibrium constant expressed in terms of activities and is a dimensionless ratio. 

From AG = AH - TAS, AG° = AH 0 - TAS° 

Substituting into AG° = -RT In K., 

and re-arranging _ A U° A S° 

n ^ —— nr + This is termed as van't Hoff ISOCHORE 


For condensed phases, Ka is expressed in terms of concentration or mole fractions, in which case AH = 
AU and, 


In K=- 


AU° 

RT 


+ 


AS° 

R 


3.1 Chemical equilibrium and metal extraction from aqueous media - 


3.1.1 Application of equilibrium constant to solubility 

The equilibrium solubility of a metal hydroxide in an aqueous acidic solution may be represented by the 
equation: 

Me(OH) 2 + 2H + = Me 2+ + 2H 2 0 (you will see later that this equation has application in 

solution purification, geochemistry, environmental 
chemistry as well) 

At the same time some of the water molecules ionize according to: 


H 2 0 = H + + OH' K w = [H ] [OH ] = 10' 14 at 25°C 


The first equation can be written as a sum of two reactions: 



with equilibrium constant K' 


1) Me(OH) 2 = Me 2+ + 20H' 

2) 2H + + 20H' = 2H 2 0 with equilibrium constant K w 

The equilibrium constant expression is: 


. = [Me^] = [ +2] 

[H +1 f L J 


[oir 


[OH ~ 1 \ [H + 


1 12 


= K 


Ki 


In general for Me(OH) z , K* = K7 K w z 
At equilibrium, 


log [Me 2+ ] e = -2pH + log K* 

The equation shows that the equilibrium concentration of the metal ions, and therefore the solubility (or 
leachability) of the metal hydroxide depends on the pH. If the concentration of the metal ion, [Me 2+ ], is less 
than the expected equilibrium value at a given pH (ie., log Q < log K* ) then the reaction proceeds to the 
right, that is dissolution occurs (or in terms of stability [Me 2+ ] is more stable than Me(OH) 2 . If on the other 
hand if metal ion concentration is more than the equilibrium value, then reaction proceeds to the left, that is 
precipitation of Me(OH) 2 occurs ( Me(OH) 2 is more stable than [Me 2+ ]). 

log [Me 2+ ] < -2pH + log K* , [Me 2+ ] is more stable. 

log [Me 2+ ] > -2pH + log K*, Me(OH) 2 is more stable. 

In graphical form, log [Me 2+ ] > -2pH + log K* is an equation of a straight line (y = mx + c), for which, 

1) m = -2 

2) at log [Me 2+ ] = 0, pH = Vi log K*. 

An example is the Fe(OH) 3 - H20 system: 

Fe(OH) 3 + 3H + = Fe 3+ + 3H 2 0, for which log K* = 3.3. 

Therefore: log [Fe 3+ ] = 3.3 - 3pH. 



The graph of this is shown 
indicating the stability regions of 
the [Fe 3+ ] and Fe(OH) 3 . 

These kinds of diagrams find 
applications in determining the 
conditions for leaching, 
precipitation, etc. 




3.1.2 Pourbaix (Eh-pH) Diagrams 

Every oxidation and reduction system can be represented by the equation: 
Re = Ox + ne 

The corresponding half-cell potential is given by: 


E h is the half-cell potential relative to the H 2 -electrode. Therefore for the half-cell reaction: 


aA + bB = cC + dD + ne , 


Ef,.25°C~ E + 


0.0592 


[cj[pY_ 

[A} a [B} b 


3.1.2.1 The stability of water 

The stability of water can be expressed as a function of hydrogen and oxygen partial pressures from the 
equilibrium constant of the equation: 


2H 2 0 ( i) - 2H 2(g) + 0 2(g) , 


giving the upper limit of water stability as the equilibrium between water and oxygen at 1 bar (atm) 
pressure. The stability as functions of Eh and/or pH can be obtained by rewriting the equilibrium equation 
between water and oxygen and/or electrons: 


2H 2 0 ( i) - 0 2 ( g ) + 4H + aq + 4e 




1 i4 


„ r o 0.0592, Po{H + } 
Eh=E 


or 


E h =E°+ 


0.0592 


log [H + 'T ( 


This simplifies to: 

Eh - E° - 0.0592 pH (The upper limit of water stability, po 2 = 1 bar, and the activity of 

pure water is 1) 

The plot of Eh - E° - 0.0592 pH is a straight line with a slope of -0.0592 Volts per pH unit, and has an 
intercept of E°. To obtain numerical values of E°, we can use the expression: 

r o A G° 

E = —E” 
nr 

Form equation (**), 

AG° r = AG°f 02 + 4 AG°f H + - 2 AG°fH 2 o 

= 0 +(4x0) - (2 x -56.69) = 113.4 kcal (....kJ) 


Therefore E° = (113.4) / 4 x 23.06 = 1.23 Volts 



Therefore Eh = 1.23 - 0.059pH 


(see Figure) 


0.059 pH 
0.059 pH 

From 2H 2 0 ( i) = 2H 2(g ) + 0 2 ( g ) 


H 200 - 

2H + aq + 2e 



In this case 

a g: 

0 



n 

a 

II 

AGT = AG° fH + - AG ° fH2 = 0-0 = 0 

Therefore 

„ 0.059 

E„ = 2 

log Ph 2 ~ 

0.059pH 


3 . 1 . 2.2 Types of Boundaries in Eh - pH space 

(1) For pH-dependent only 

M 2 0 3 + 6 H + = 2M 3+ + 3H 2 0 p H - - i og 

6 

(2) For Eh-dependent only: metal to ion 

M =M + + e- Eh = E° + ^^logM +1 

n 

(3) For Eh-dependent only: ion to ion 

Eh = E ° + = E° 

M 2+ = M 3+ + e" 77 ] 


In general equation (***) can be written as: 
or in terms of hydrogen partial pressure: 


r 0.059. 

E h = 1-23 + log p 02 - 

r r° °- 059 1 

E h = E -^log^, - 


(4) For Eh-pH dependent: activities of redox species equal 
M 3+ + 2H 2 0 = MOT + 4H + + 2e‘ 


„ . 0.059, [MO +v 

E h = E + -log^—— 

n [M ] 


+ ^^ log [// +1 ] 4 


E h = E°- 0.118 pH 

(5) For Eh-pH dependent: solid - ion redox reactions 
MO + 2H + = M 3+ + H 2 0 + e" P 

h 


0.059, rj/r+3l 
+ -log [AT J ] 


0.059 


n 


log 


[H 


+ 1 12 



E° + ^y^log[M +3 ] + 0.0295 pH 


The slope of this boundary is positive since hydrogen ions 
appear with reactants. 


(6) For Eh-pH dependent: solid - ion redox reaction 

MO + h 2 0 = M0 2 + + 2H + +e- E h = E° + log[ MOt 3 ] + 0.118 pH 

The slope of this boundary is negative since hydrogen ions 
appear with the products. 

(7) For Eh-pH dependent: solid - ion redox reaction 


MO + H 2 0 = M0 2 + + 2 FT + e 


E h = E° + ^piog[MO_: 3 ] + 0.118 pH 


(8) For Eh-pH dependent: more than one element involved in redox reactions: 

a) Oxidation of metal sulphide: 

2MS + 11H 2 0 = M 2 0 3 + 2S0 4 2 ' +22H + + 18e' 


E h = E° + 


E, = E° + 


0.059 

n 

0.059 


log [SO. 


-2i2 


0.059 


log[N + 


1 122 


18 


log [ SO 4 2 ] + 0.072 pH 


b) Oxidation of M ion and S(-II) sulphur 


2M 2+ + 2H 2 S + 11H 2 0 = M,0 3 + 2SOT +26H + + 18e 


E„ = E° + 


0.059, [SO~ 4 -]- 0.059, r „ + M26 

lo §T—ITT + - l °gi H J 


n 


'M 


n 


1-2 1 


Z7= A 0.059 , L so ; j 

Eh = E + —z— log-— -j-j - 0.085pH 


18 


[M + 


c) S oxidized, M reduced 


2MS + 4H 2 0 = M° + S0 4 2 ' +8H + + 6e‘ 


E h = E° + M59 log[iSC >- 2 f 


M^iogt ^ 1 ] 8 



E, — E° + ^^log[S<9 4 2 ] + 0.079 pH 
6 


3.1.2.3 Preparation of Stability Diagrams (Eh - pH, [M ] - pH, [L] - pH) 

The stability diagrams are made up of stability lines and predominant areas. To prepare stability diagrams 
you need: 

1. Thermodynamic data 

2. Mathematical manipulations 

Step 1: Write the relative stability equations using system variables (e", H + , [M], H 2 0). For example for 
M 2 0/M 2+ system, the following equation is used: 

2 M 2 

[H +1 f 

log Q = 2 log [M 2+ ] + 2 log (e ) - 2 log [H + ] 

= 2 log [M 2+ ] - 2 p£ + 2pH 


M 2 0 (s) + 2H + = 2M 2+ + H 2 0 + 2e" , for which 


Q=[M +1 


Step 2: Get the equation for the relative stability lines 
Step 3: Plot the relative stability lines 

- the reaction goes to the right if logK > log Q; (log K > 2 log [M 2+ ] - 2 p, + 2pH) 


- at equilibrium, log K = log Q (log K = 2 log [M 2+ ] - 2 p£ + 2pH) 


P 


pH 

Step 4: Identify the relative stability regions 

The relative stability of M 2 0 or M 2+ is determined by three variables; {M 2+ }, (e), (H + ). There are 
three possibilities: 

(a) for {M 2+ } is constant, then: p£ = {-Vi log K + Vi log [M 2+ ] }+ pH 

We know that the reaction proceeds to the right if: 



log Q < log K or on rearranging, 


p> {- l A log K + !/£ log [M 2 +]} + pH 




That is if p£ > k + pH, ( k = const.) then M 2+ lies above the equilibrium line 

(b) for p£> k the reaction proceeds to the right if: 

log K > 2 log [M 2+ ] - 2p£ + 2pH or on rearranging, 

log [M 2+ ] < [ Vi log K + pe] - pH or log [M 2+ ] < k - pH 

M 2 0 lies above the stability line 



(c) for pH = constant, p£ > k + log [M 2+ ] the diagram is similar to that of (a) 


3.2 Chemical equilibrium and the stability of metal oxides in aqueous media 

The stability relations among various metal oxides of metallurgical interest can be obtained by writing 
reactions in terms of the metal compounds and gaseous oxygen, including all permutations and 
combinations possible. Calculating the equilibrium constants and knowing that the activities of crystalline 
solids at standard conditions is 1, the equilibrium constant is written as a function of oxygen pressure. 
Calculating the standard free energy of each reaction and using, 

AG° = - RT InK 

the partial 02 pressure can be determined, and from this the relative stabilities of the oxides. 

In the case of iron the following reactions could be taken into consideration: 

Fe (s) + ViOzfg, = FeO (s) (wustite) 


3Fe (s) + 20 2 ( g) = Fe 3 0 4 (s) (magnetite) 





(hematite) 


2Fe (S ) + 3/20 2 ( g ) - Fe 2 0 3 ( S ) 

3FeO (s) + V40 2(g ) = Fe 3 0 4 (s) (magnetite) 

2FeO (s) + 1 /i0 2(g) = Fe 2 0 3(s) (haematite) 

2Fe 3 0 4(s) + ViOn^ = 3Fe 2 0 3 ( S ) 

From the above, it could be seen that the most stable iron phases are native iron, magnetite and haematite. 
In aqueous media, the stability relations can be converted to relations expressed as functions of E h and pH 
by adding water dissociation half-cell to the reactions, eg. oxidation of Fe to magnetite: 

3Fe (S ) + 20 2(g) — Fe 3 0 4 ( S ) 

4H 2 0(1) = 8 H + aq + 20 2 + 8 e 

3Fe (s) + 4H 2 0(1) = Fe 3 0 4(s) + 8 H + aq + 8 e 

A G° 

nF gives E° = -0.084 V 

or 

E h = -0.084 + M^log[// +1 ] 8 ; Eh = -0.084-0.059pH 
8 

Similarly the oxidation of magnetite to haematite can be represented by: 

Eh = 0.221 - 0.059pH 

The same could be done for iron hydroxides, ions in equilibrium with iron oxides, etc., and a composite 
diagram could be composed showing the stability regions of the various species in the Eh-pH diagram. 

(GIVE AN EXAMPLE OF THE Cu - H 2 0 SYSTEM) 


E° 


Calculating the )G° r and substituting into 


„ n n « 4 j. 0.0592. [Fe,0,][H"} 

E h - -0.084 + — 7 —log 


8 [ Fef[H,0 ] 4 


3.3 Chemical equilibrium and the stability of metal oxides in pyrometallurgy 

The relationship between the entholpy change, AH (heat of reaction), entropy change, AS, and the free 
energy change AG. 

AG = AH - TAS. 

At standard state (pure condensed phase or 1 atm. pressure for gases): 



AG° = AH 0 -TA° S. Sometimes AG° is given as a function of temperature, T, in which case 
AG° = AH 0 + a T log T - % abT 2 - 1/6 cT 3 + IT (a, b, c , I are constants) 

For reactions at equilibrium: 

G (T, p) = O, and 

AG° = RT In K K: equilibrium constant (so AG enables the calculation of K and vice versa) 

For a metal-metal oxide system: 

2 

r\ r\ 

— Me + 0 2 ——Me x O and A G- - RT In K = - RT hi 

y y 

a Me 

(Note that Me and Me x O y are assumed to be in standard state) 

Combining with the above equation gives: RT In p 0 A G° — A // ' — T AS" 

We can therefore say that with good approximation, AG° or RTlnPo^ has a linear function with 

temperature. Plots of AG° versus T are known as Ellingham diagrams. It is also possible to regard the 
metal-metal oxides as being in equilibrium with a gas mixture, which in itself defines the oxygen partial 
pressure (eg. a mixture of CO, CO 2 ) 

Thus for: 

Me + C0 2 = MeO + CO AGj 0 = AHj° - TASj 0 = - RT In (pC0/pC0 2 ) 

This is a difference of two reactions: 

(1) 2Me +0 2 = 2MeO AG 2 ° = AH 2 ° - T AS 2 ° = RT In p0 2 

(2) 2 CO + 0 2 = 2C0 2 AG 3 = A H\ - TRS\ = RT\n{ pC0 ^ l pC0 ) 2 

p0 2 

so that, 4G° = [AH° 2 - AH°) - T{AS 0 2 - AS° 3 ) 

It is possible by these to plot AG = f (T) for various oxides as presentd in the Ellingham diagrams. Detailed 
discussion of the Ellingham diagrams will be don under pyrometallurgy. 





































































































































































4 COMPLEX EQUILIBRIA 

Another application of chemical equilibrium is in the formaion of complexes Metal ions in aqueous media 
do not exist as bare ions, but are present as complexes (coordinated compounds). In the absence of complex 
forming ligands, the central metal ions are surrounded with water molecules through coordination forming 
a kind of loose bond with the ion. These "hydrated" ions are called aquo-complexes. When other complex 
forming ligands are present in the solution then depending on the thermodynamics of the complex 
formation (or the stability of the complexes) the water molecules are replaced by the complex forming 
ligands. 

The complex forming ligand could either be an ion (anionic ligand) or a molecule (neutral ligand) and the 
number of ligands that could be bonded to the central metal ion (nucleus) depends on the coordination 
number of the nucleus. Discussions on coordination numbers can be found in text books of inorganic 
chemistry. The ligands are bonded to the central metal ion by either ionic or covalent bonding. 

If you consider for example the complex, [Co(NH 3 ) 6 ] 3+ .3CF, the ammonia molecules are needed to satisfy 
the coordination number (6) of the central metal ion, whilst the CT ions are needed to satisfy the oxidation 
number (3) of the central metal ion. 


4.1 Complexes with Neutral Ligands 

Neutral ligands are those molecules with single donor atoms, which have no charges, eg. H 2 0, NH 3 . An 
example is the formation of Cr 3+ amine complex: 

[Cr(H 2 0) 6 ] 3+ + NH 3 = [Cr(H 2 0) 5 NH 3 ] 3+ + H 2 0 

[Cr(H 2 0) 5 NH 3 ] 3+ + NH 3 = [Cr(H 2 0) 4 (NH 3 ) 2 ] 3+ + H 2 0 

[Cr(H 2 0) 4 (NH 3 ) 2 ] 3+ + NH 3 = [Cr(H 2 0) 3 (NH 3 ) 3 ] 3+ + H 2 0 

[Cr(H 2 0) 3 (NH 3 ) 3 ] 3+ + NH 3 = [Cr(H 2 0) 2 (NH 3 ) 4 ] 3+ + H 2 Q 


[Cr(H 2 0)(NH 3 ) 5 ] 3+ + NH 3 = [Cr(NH 3 ) 6 ] 3+ + H 2 0 
The overall complex formation reaction is: 

[Cr(H 2 0) 6 ] 3+ + 6NH 3 = [Cr(NH 3 ) 6 ] 3+ + 6H 2 0 

The above represent substitution reactions (substitution of water molecules with ammonia molecules), 
taking place in steps, with the formation of intermediate mix ligand complexes. For the amine complexes to 
be formed the concentration of ammonia in the solution must be high enough. In this case other side 
reactions may also take place, notably: 

NH 3 + H 2 0 = NH 4 + + OH- 

As a result, the formation of insoluble hydroxide precipitates may also occur, eg. 

Fe 3+ + 30H- = Fe(OH) 3 


However, if the hydroxide is soluble and the amine complex stability is high then the formation of the 



amine complex is more probable. For example, in the case of Cu 2+ , Cu(OH) 2 is soluble and the reaction in 
presence of NH 3 therefore is: 

Cu 2+ + 4NH 3 = [Cu(NH 3 ) 4 ] 2+ 


4.2 Complex formation with anionic ligands 

- acido complexes: the ligands in this case (fluoro-, chloro- bromo-, and iodo-complexes) are 
slightly basic and therefore their formation depends slightly on the pH. 

- ciano complexes: are formed only in neutral or basic medium, with the cyanide ion replacing the 
water molecules around the metal ion, eg. 

[Fe(H 2 0) 6 ] 2+ + CN- = [Fe(H 2 0) 5 CN] + + H 2 0 


[Fe(H 2 0)(CN) 5 ] 3 " + CN" = [Fe(CN) s ] 4 - + H 2 0 

Other cianocoplexes are [Fe(CN) s ] 3 ~, [Au(CN) 2 ]", [Ag(CN) 2 ]', [Zn(CN) 4 ] 2 ', [Ni(CN) 4 ] 2 '. Some of 
these complexes can also be formed in weakly acidic solutions since their decomposition is 
kinetically inhibited. 

-hydroxo complexes are formed from aquo-complexes by the deprotonation of the water molecules 
around the metal ion in the presence of a proton acceptor, ie: 

Me(H 2 0) m+ n = Me(0H)(H 2 0) (ml)+ n .i + H + 

The role of the proton acceptor can be played by the water molecule since the coordinated water 
molecule is a stronger acid than the non-coordinated. If the pH of the solution is increased by the 
addition of NaOH, then all the coordinated water molecules are exchanged for OH" ions: 

[Me(H 2 0) n ] m+ + H 2 0 = [Me(H 2 0) n .i(OH)] (m ' 1)+ + H 3 0 + 

[Me(H 2 0) n . 1 (0H)] (ml)+ + H 2 0 = [Me(H 2 0) n _ 2 (0H) 2 ] (m - 2)+ + H 3 0 + 


[Me(H 2 0) n . (n _ 1) (0H) (n .i ) ] m ‘ (n ' 1)+ + H 2 0 = Me(OH) n + H 3 0 + 


Even though in all cases the initial and final results are known, the intermediate steps are not often 
the same as above. Polynuclei complexes may be formed, especially in the case of amphoteric 
hydroxides, eg.: 

2A1(H 2 0) 3 (0H) 3 = A1,0,.9H,0 


Or in the case of Cr(III), the following represents the complexation route: 



Cr{H 2 0\(0H\ 


+ 3¥L 2 0 


+ 3H 2 0 


[Cr(H 2 0 \] +3 + 30H 1 [Cr{H 2 0\] +i + 3H 2 0 + 


Other complexes include: 

-chelate complexes, eg. with EDTA (ethylene diamine tetraacetic acid) 
-polynuclei complexes, ie with more than one central atom taking part, 
-mix ligand complexes, eg [PdCl 4 ] 2 ", [PdCl 3 Br] 2 ", [PdCl,Br,] 2 ‘ 


4.3 Complex Stability Constants 

The thermodynamics of complex formation is shown through the relationship between the equilibrium 
complex formation constant (K) and the free energy change during the process. 

Considering the equilibrium complex formation reaction between metal ion (Me) and ligand (L): 

Me + L = MeL, for which from mass action law, 

a MeL K is the consecutive stability constant. 

K= - a, 

a L 

u Me 

If the reaction proceeds through intermediate complexes, then: 


Me + L = MeL 


K x = 


a 


MeL 


a Me a L 


= fix 


Me + 2L = MeL, (MeL + L = MeL 2 ) 


K 2 = 


MeL, 


a MeL a L 


; P2 


MeL- 


a Me a L 


2 = K x . K 2 


Me + 3L = MeL 3 (MeL, + L = MeL 3 ) 


K 3 = 


1 MeL, 


a MeL, a L 


; fit = 


1 MeL, 


a Me a L 


= K 1 . K 2 . K 3 


Me + nL = ML n (MeL n _i + L = MeL n ) 


* * 
fin = n K, 


(3 is known as the overall stability constant. The physical understanding of the consecutive stability 
constant is that, in a given system containing metal ion (Me) and ligand (L), all the species from MeL ... to 
MeL n exist and the amount of each species depends on the concentration of Me and L, and on the individual 
stability constants (K,). For example: 












Cd 2+ + NH 3 = [Cd(NH 3 )] 2+ 

Ki = 10 2 ' 65 

[Cd(NH 3 )] 2+ + NH 3 = [Cd(NH 3 ) 2 ] 2+ 

K 2 = lO 210 

[Cd(NH 3 ) 2 ] 2+ + NH 3 = [Cd(NH 3 ) 3 ] 2+ 

K 3 = 10 144 

[Cd(NH 3 ) 3 ] 2+ + NH 3 - [Cd(NH 3 ) 4 ] 2+ 

k 4 = io°- 93 , p 4 = 10 712 

the case of cyanide complex, 

Cd 2+ + CN' = [Cd(CN)] + 

Ki = 10 548 

[Cd(CN)] + + CN' = [Cd(CN) 2 ] 

K 2 = 10 s 12 

[Cd(CN) 2 ] + CN’ = [Cd(CN) 3 ]" 

K 3 = 10 463 

[Cd(CN) 3 ] + CN' = [Cd(CN) 4 ] 2 " 

K 4 = 10 355 , (3 4 = 10 188 


The above shows that in such a system MeL is first formed, and as the ligand concentration increases MeL 2 
increases whilst MeL decreases, and so on. 

Thermodynamic Tables exist where the overall stability constants are listed for various number of 
complexes. Sometimes the logK is listed. 

The stability of complexes for a particular metal ion depends on the size of the ion, the charge on the ion 
and some steric factors. Generally, complex stability constant decreases in the following order: 

Li + > Na + >K + >Rb + >Cs + ; 

Mg 2+ >Ca 2+ > Sr 2+ > Ba 2+ ...; 

Al 3+ >Sc 3+ ... 

Mn 2+ < Fe 2+ < Co 2+ < Ni 2+ < Cu 2+ < Zn 2+ 


4.4 Complex formation in the presence of other reactions 

It is possible that complex formation takes place together with other reactions. Such reactions may be acid- 
base, redox, precipitation or other complex formation reactions. The diagram below shows reactions that 
are possible. In such cases the nature of the reactions determine whether the complex formation reaction is 
the main or side reaction. 
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4.5 Complex Formation and Solubility 

In some applications in the extractive metallurgy the formation of complexes is beneficial because it leads 
to the stability of the metal ion in solution and therefore enhances the solubility of the metal during 
leaching. In other applications the stability of the metal complexes could be used in selective separation of 
different metals from one another. 


A typical example is the leaching of gold in cyanide solution. As native gold is oxidized in the presence of 
oxygen, the aurous ions forms a complex with the cyanide ions present in the solution: 

Au = Au + + e E° = 1.69 V 

Au + + 2CN’ = [Au(CN) 2 ]- (3 = 10 38 (mol/dm 3 )- 2 


E 


r0 RT . [Au[CN) 2 ] 1 

E + -In--- 

nF p[CN\ 2 


Combining all constants with EE and rearranging, the standard electrode potential for complex formation 
(E° C omp) is obtained: 


E 


O 

comp 


E° 


RT , n 

rF n P 


1.69 - 2.24 = —0.55V 


Therefore, 


E 


E° 

comp 



nF 


AulCN 

y 

.CAT 7 

2 


+ 




The above shows that because of complex formation the activity of the metal ion decreases leading to a 
shift in the electrode potential. 


Another example is the change in solubility as a result of the addition of complex forming ligands into the 
solution: 

Cu 2 S + 8CN’ = 2 [Cu(CN) 4 ] 3 " + S 2 ' 

This particular reaction is undesirable in gold cyanide leach systems because it consumes cyanide. 

4.6 Speciation Diagrams 

Using data from the discussions under the section on “Complex Stability Constants”, it is possible to draw 
speciation diagrams indicating the regions of stability of the complex species. 

Let the total concentration of dissolved metal be C M 

Then C M = [Me] + [MeL] + ... + [MeLJ 

Cm = [Me] + PJLUMe] + P, [L] ! [Me] + Pi [L]' [Me] 


C M = [Me\ { 1 + £ /?, [if) 

i=l 




i= 1 


The speciation diagram can then be drawn by plotting: 

(i) a! vs log [L] 

(ii) 0U vs pL (- log [L]) 


or a f = [L] 1 a Q 


Example : the speciation diagram for Ag - OH - H 2 0 system 
Ag + + H 2 0 = Ag(OH) aq + H + 

(log $, = -12.0) 

Ag + + 2H 2 0 = Ag(OH) 2aq - + 2H + 

(log $ 2 = -24.0) 




5 PHASE EQUILIBRIA 

5.1 Introduction 

The knowledge of phase equilibria is important not only in the study of phase processes taking place at 
ordinary atmospheric conditions, but also in the study of high temperature reactions, processes involving 
different phases. These include high temperature and pressure solid-solid, solid- liquid, solid-gas and 
liquid-gas reactions, the development of products in the ceramics and glass industries, alloys and special 
alloys, the formation or precipitation of minerals phases, among others. 

5.2 The Thermodynamic Phase Equation 

The complexity in the study of such systems increases as the number of components and phases increase, ft 
is easier in this case to first of all understand simpler multiphase, multicomponent systems and then apply 
the principles to the more complicated systems. 

Consider a one-component system involving two phases (ot , (3). 

The two phases are in equilibrium if the following conditions are fulfilled: 

Ta — Tpj p« - pp; and Pa — M3 

From definition, p« = g«, where g- is the molar free energy. 

Similarly M 3 - 93 

At equilibrium and from the expression, dG = -SdT + Vdp, 

dga = -SadT + Vadp = dgp = -spdT + Vpdp 

or, (v a - v^)dp — (Sa-sp)dT 





dp _ (Sq-Sp) 

dT (v a Vp) 


For isothermal transformations and at equilibrium, 

Ag=0=Ah — T As 


This is known as the Clausius-Clayperon Equation. This equation 
shows the relationship between pressure and temperature at which 
phases " and IS exist in equilibrium. It also shows the enthalpy and 
volume changes that accompany such changes. 

For liquid/vapour considerations the volume of a unit mass of the condensed phase is negligible as 
compared to that of the vapour phase, and assuming ideal conditions therefore, 


AG=0 = AH - TAS 



Therefore: dp _ Ah 

~dT~ T Av 


Av- 


RT 


P 


dp 

~dT 


In the form of a differential equation, 


p Ah 
RT 2 


dlnp 

Aj) 


Ah 

~R 


Given points such as the melting and boiling points in the system, the equations can be integrated and 
presented as a plot of p against T, or T against % composition. 


5.3 The Phase Rule 

For multicomponent multiphase systems the following conditions must be fulfilled for equilibrium to be 
established among the phases: 

Tq = Tp = T y = ... (P-1 equalities, P: the no. phases) 

p« = pp = p Y = ... (P -1 equalities) 

Micx= MiP = Miy = ••• (P-1 equalities) 

M23 = M215 = |J2y = ... (P -1 equalities) 


In this type of system, there are C+l variables per phase, and therefore P(C+1) variables in P phases. At the 
same time there are (P-1) equalities in chemical potential for each C component, satisfying (P-l)(C+2) 
equations. 

Therefore the Degree of Freedom (F) is given by: 



F = P(C+1) - (P-l)(C+2) 

(variables) (equalities/equations) 

F = C — P + 2 This is the GIBBS Phase Rule 


5.4 Phase Diagrams 

In a PVT system, the intensive state is established by fixing the temperature, pressure and composition of 
the phases in the system. In equilibrium these variables are not independent, and by fixing a limited number 
of them the others are automatically established. For example for a mono-component system like water, 
once the pressure is fixed the boiling or freezing point is automatically known. 


Examples of Phase Diagrams: 
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6 MIXTURES, SOLUTIONS AND MELTS 


6.1 General Principles 

Substances often form a variety of mixtures. The components (metals, elements, compounds) often keep 
their individual properties. These mixtures can be homogeneous (eg. solutions) or heterogeneous. It is 
therefore necessary to study the relationship between the compositions structures and properties of these 
mixtures. 

When a system is made up of a mixture of different species (eg. chemical species), one would expect the 
mixture property (M, thermodynamic property, eg. U, S, H, etc.) to be related to the properties (M) of the 
pure species which make up the mixtures. There are, however, no generally valid relationships. In any 
particular system the connection must be experimentally determined. 

In a mixture or solution the constituents do not have any identifiable separate thermodynamic properties. 
However, it is good to think of the mixture property as being apportioned among the mixture constituents in 
some way. Once this appropriate procedure is adopted, the resulting property values may be treated as 
though they were properties of the individual constituents as they exist in the solution. 

The total moles of mixture is given by: 

Then the total system property (nM) is: 

where M is the partial molar property. It denotes how the mixture property is apportioned among the 
mixture constituents. If the resulting molar properties are treated as properties of constituents in solution 
then it can be said that the total mixture property is a sum of the properties attributed to the constituents. 

nM = / M = X(-M) = XU^) 

n 


« = X n i 

nM - Xk^) 


Differentiating the above expression gives: d{nM) = X («, d M ,) + X M, d n, 


The molar properties can be expressed as a function of measurable quantities such as temperature, pressure 
and the number of moles in the system: 

nM = f{T,p,n t ) 


which on differentiation gives: d (nM 


d ( nM) 
dT 


dT + { d ^h T .„ d P + X M<dn t 


The expressions for d (nM) are true only if: 
/ d ( nM ), ,d{nM 


QJ, Jp.n dT + V QJ, )T 


_ This is the general form of the Gibbs- 

,n dp + 2-iMjdnj — O Duhem equation 


or in terms of mole fractions: 


d(nM , 
dT 


dT + 


d ( nM ) • 
dT 


dp + X M i dx i —0 


An application of the Gibbs-Duhem equation is the contribution of the various species to the overall molar 
free energy (chemical potential) on mixing: 



dT + (^)r,, dp + £*,</G,=0 
In other words: _ sdT + Vdp~Yj x t d\i t =0 

which at constant T and p simplifies to: I Xjd m=o 

The quantitative study of complex real solutions can be provided through the study of ideal gas mixtures. 
For an ideal gas: 

pV = nRT where n is the number of moles of gas and R is the gas constant (R = 8.314 
J/mol.K). 

For a mole of ideal gas in a closed system and at constant temperature: 
dG = Vdp 

dG = (RT/p)dp or dG = RT dlnp, for which 

G 2 - Gi — RT ln(p 2 /pi), giving only the change in free energy. 

Generally, 

G(p, T ) = G°(d + RT lnp , where G 0 m is the free energy at standard state. 

Similarly, 

M(P,T) = M°m + RT lnp 


6.2 The change in system property on mixing 

Before mixing the thermodynamic property of the components is given by: M t )=^ x t M 

The mixture property is expressed as: ,y/ ( X; Mj ) 

The change in thermodynamic property on mixing can therefore be calculated as: 


A M . -M, 


-M, 


pure 


&M mL =Y; XiM-Yj x i M 'i 




For ideal gas mixtures. 

PiV = n,RT (Dalton's Law) 

Pi = XjP, where P is the total pressure (P = X p, and n = Zn,) 

Mup.t) ~ M i(T) + RT lnp; = p i(T) + RT lnx,P 
On mixing, 

G (m«) = S = Z n^. + RT X n M Pi 



Before mixing the free energy of the pure component is: 

The free energy change on mixing is: A G {mix) - X n,(ii, 


G( pU recomp)=ll n^ + RT X Pi 

-[ii)=RT^nil ny 

1 i 


A special case is if volume is constant (V = const.), in this case, 

Pi = P 2 = P 3 = ... - P, 

Then: A G [mix) = RT^ n t In x, 

For a mole of gas the above equation simplifies to: ^ — pp X x In x 

Similarly the entropy change on mixing is: A H A G 

J r j o o AC _ mix mix 

^ ^ (mix ) rji rj-f 

Simplifying, \S mix = -R^x^nx, 


For real gas mixtures 

The general gas equation is not valid. In this case the equations that may be applied are the van der Waals 
and Varial equations, as well as the equation containing the compressibility factor (pV = zRT). However, to 
maintain simplicity we introduce a new thermodynamic constant called fugacitv (f) , which is defined as: 

lim — = 1 

/>-»o p 

So for a real gas therefore, p (PiT) = p" (T) + RT In f (p , T) 

And for real gas mixtures, p i( p, T) = p° im + RT In f i( p >T ) 

6.3 Aqueous Solutions 

Leading from the discussions on gases, aqueous solutions can either be ideal or real. A simplified definition 
of an ideal solution is one in which the cohesive forces in the solution would be the same as those existing 
in the separate components of the solution. 

When two substances (A & B), whose molecules are similar, form an aqueous solution the vapour pressure 
of the mixture is related to the vapour pressure of the pure components by the relation: 

Pa = X A Pa* and p B = x B P B *, where 

p is the equilibrium partial pressure, P is the vapour pressure of pure substance and x is the mole fraction. 
The above relation is known as Raoult's Law, which is generally stated as: 

Pi = X; P;* Raoult's Law 

Solutions that obey Raoult's Law are said to be ideal solutions. When the Law is applied to the solubility of 
gases in liquid, then it becomes known as Henry's Law, which is defined as the concentration of any single 
molecular species in two phases in equilibrium bear a constant ratio to each other at a fixed temperature: 


- 1 

Xi-yPi-kPi 



Henry's Law 


Hx i -p i 


The above equation shows that the mole fraction of a gaseous solute is proportional to the partial 
pressure(pi) of the gas above the solution. As P* increases with temperature, the value of k decreases. 

Raoult's Law is a special case of Henry's Law and it is only valid for dilute solutions. A limiting case is 
represented by pure solvents (molecular interactions are the same) or very similar solvent and solute cases. 


(Diagrams) 



6.4 Representation of Solution Concentrations 

The composition of solutions are usually expressed in terms of: 

1) molarity (M) - moles of substance per litre of solute 

2) molality (m) - moles of substance per 1000 gram solute 

3) mole fraction (x) - ratio of moles of a component to the total number of moles in solution. 


For binary solution: 

n i n 2 

x,=- x 2 =- 

n x +n 2 n l +n 2 


For binary aqueous solutions: 


For an M; molar solution: 


n : 


n,+n 


H~0 


M. 

l 

(lOOOp —M,TF,) 
M * +1 — 


because the number of moles of water in a litre of 
solution is (1000D - MiW;), where W; is the molar mass 
of component i. 


For infinite dilution the relationship between molarity or molality and mole fraction is given by: 
18x. 18 x. 

M = - ' 77 ? =- 

' 1000p ’ ' 1000 


Another standard state often used in the case of liquid (molten) metals is the 1 weight percent standard 
state fl wt%). This is located on the Henry's Law line at 1 wt%. 


6.5 Ideal and Non-ideal Solutions 

Thermodynamically, at equilibrium the chemical potentials of the components of solutions that follow 
Henry's Law is given by: 


Mi (Uq) Mi (vap) M i * RT In Pi 


Mi (Uq) = M°i + RT In ki + RT In X; = M*i + RT In X; 
p*i is the standard state chemical potential in the liquid phase. 

M°i is the standard state chemical potential in the vapour phase 

Similar to the discussions involving real gases, for real solutions we introduce another thermodynamic 
property called the activity (a) . By definition the activity (aO of a component in any state at temperature, 
T, is the ratio of the fugacity of the substance in that state to its fugacity in its standard state: 



Another probably more practical definition is the "effective concentration of ions, atoms and molecules". 
This is introduced so that in non-ideal cases it will be possible to calculate. Activity is calculated as the 
product of the mole fraction (x) or concentration (c) and what is known as the activity coefficient (y): 
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a i =x i y i 

For real solutions therefore: 


Hi = |j*i + RT In a; = [j*; + RT In X; + RT In Y; 
At equilibrium, 

Mi(vap) = H' i + RT In r = Miciiq) = M*i + RT In a; 


a i~Pi exp 


Uy-ih) 

RT 


On mixing ideal gases and ideal solutions, the enthalpy and volume changes are zero. The free energy 
and entropy changes on mixing ideal solutions are: 

A G=X n iP~ X n iR( — X n i RTlnx i 

In terms of mole fractions, AG= RT^ x \nx 


The entropy change on mixing is given by: 


AS mi =- 


A H 


A G„ 


- = — i? ^ x, In x f 


Unlike in the case of ideal solutions, mixing real solutions is accompanied by enthalpy and volume 
changes. In addition the free energy change is given by: 

A G mix =RT X x,ln a t =RT X x, ln x t +RT X x i \ny i 


The first term represents the change due to ideal solutions and the second term represents excess free 
energy per mole of solution on mixing. 

In general, G = G id + G xs and a n —\G id +G XS 

f-j l—S. VJ • VJ • VJ 


For any solution, A G mix = A H mix - T A S imx 

G xs = A G mix -AG id mix =AH mix -T{A S mix -A S id mix ) 

For binary solutions, G xs = RT[x A \ny A +x B \n y B ) 


The diagrams showing the change in free energy, entropy, enthalpy and volume are given below: 


diagrams 
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6.6 The Tangent Intercept Method 

This method enables us to determine the partial molar quantities of the various thermodynamic properties 
(G, H, S, V, etc.,) of the components making up the solution. 



6.7 Solutions of Electrolytes 


6.7.1 Basic Concepts 

Electrolytic solutions are special form of solutions in which the solute exists partially or completely in the 
form of ions. By definition an electrolyte is any compound which in solution or fused state conducts 
electric current and is decomposed by it. 

The behaviour of electrolytes has been studied by notable scientists including Baron C.J. von Gotthuss 
( who suggested that molecules of electrolytes were polar with positive and negative ends, and in an 
applied field line up in a chain with the molecules at the end dissociating and the free ions formed 
discharged at the electrodes), Classius ( who proposed that energetic collisions between undissociated 
molecules in electrolytes maintain at equilibrium a small number of charged particles and that the 
particles are responsible for observed conductivity). Kohlhausch investigated extensively conductivity of 
electrolytic solutions. Svante A. Arrhenius compiled data available at the time and proposed the theory 
for the behaviour of electrolytes. Some flaws were later discovered in his proposals especially in 
connection with the behaviour of strong electrolytes. 


P. Debye and E. Hiickel provided the theory for the modern treatment of electrolytes, that is, strong 
electrolytes are completely dissociated into ions and that the observed deviation from ideal behaviour was 
due entirely to the electrical interactions of the ions in solution. The deviations are greater for more 
highly charged ions and in more concentrated solutions.A revised theory can be summarised as: 

1) a solution of electrolytes contains free ions 

2) passage of a current through a solution of an electrolyte depends on (i) the number and (ii) the 
speed of the ions present. 

3) in a solution of a weak electrolyte the degree of ionisation increases with dilution until at 
infinite dilution there is complete ionisation. 

4) in a solution of a strong electrolyte there is always complete ionisation but ionic interference 
limits the movement of the ions. This interference decreases with dilution. 



+ 

_ + 

+ 


Ions in solution usually have an atmosphere of opposing charges 
around it. Under an electric field the ions move to the opposite 
electrodes. However, the movement is hindered by: 

1) the net accumulation of opposite charge behind the ion due to 
the asymmetric atmosphere created around the ion and exerting 
an electrostatic drag 

2) frictional drag of the solvent according to Stake's Law 

3) the ions comprising the atmosphere around the central ion 
which themselves try to move in the opposite direction under the 
applied field (Electrophoretic effect). 


According to Ohm's Law the current passing through an electrical conductor is given by: 


The resistance of an electronic conductor is given by: 
Electrical conductivity is expressed as: ^ = k 


R= p( 


/ 


v_ 

R 


In the case of solutions, electrical conduction is given in terms of equivalent conductivity (A), which is 
expressed as: 



where c' is concentration, eqv./cm 3 
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6.7.2 Migration of Ions 

Hittorf prposed that the the amont of ions in the vicinity of an electrode changes during electrolysis and 
that both anions and cations migrate in the solutions to the electrodes at different speeds carrying a 
fraction of the current through the eectrolyte. He further suggested that the fraction of current carried 
depends on their velocities. The fraction of total current carried by the ions species was called the 

Transport number (t): 

Current carried by anion Current carried by cation . 

--- —t • --- —t ■ t +t — 1 

Total current ° + ’ Total current 


From Kohlrausch's Law 


to + X 0+ and t 0 _ X 0 - 

Kohlrausch's law of independent migration (mobilities) of ions states that the quivalent conductivity 
at infinite dilution of an electrolyte (X 0 ) is equal to the sum of the ionic mobilities of the ions produced by 
the electrolyte: 


A c —X 0+ +X 0 . where X is the equivalent ionic conductivity. 


According to Arrhenius, the degree of ionization (ot) is given by the ratio of the equivalent 
conductivities at a given concentration to that at infinite dilution, 


Applying mass action principle to dilution, Ostwald obtained the Dilution Law, that is, for an electrolyte 
(AB) dissociating in solution with a concentration of c at equilibrium: 


AB = A + + B‘ 

c(l - (X) CXC 0(C 


for which the dissociation constant(K) is expressed as: 


K-- 


(a 2 c) 
1 — on 


The mobility (in m 2 /sec.V) of an ion is defined as its velocity in an electric field. For an electrolyte or 
salt (A a Y y ) which dissociates on solution in water: 


A a Y y = aA z+ + yY z ", 

the activity(a 2 ) of the electrolyte solute is found to obey Henry's Law and the value is given by: 

a 2 = k • m 2 2 for uni-univalent electrolyte solutes and symmetrical salts in 
which anion and cation have the same valency, eg. NaCl, LiF, MgS0 4 , CaC0 3 , 
etc. 

a 2 = k • m 2 3 for ternary salts such as BaCl 2 , dissolving into three moles of 
ions. 

a 2 = k • m 2 v+y for A v Y y = vA z+ + yY z ‘. 


The value of k depends on the standard state chosen. For electrolyte solute the hypothetical state is 
established by the behaviour of the solution at infinite dilution, 


a 2 = m 2 2 


as m 2 approaches 0 










At infinite dilution interaction between anions and cations vanish, and for symetrical electrolytes, 


m 2 = m+ = m_ and therefore a 2 = (m+ )(m_). 

The individual ionic activities are: 

a+ = m+ as m 2 approaches 0 

a. = m. as m 2 approaches 0 

Therefore, a 2 = (a+ )(a_ ) since m+ = m. 

These approach 1 at 
infinite dilution. 

i 

a 2 

— y± approaches 1 at 
m 2 infinite dilution. 

For unsymmetrical electrolytes, (A v Y y ) eg BaCl 2 ; 

a 2 = km 2 v+y 

a 2 = (a+ ) v (a. ) y 

a ± = [(a+ ) v (a. ) y ] 1/(v+y) 
m± =m 2 [(v) y (y) y ] 1/<v+y » 

Y± = a ± / m ± = [(y+ ) v (Y ) y l 1/(v+y) 


The mean ionic activity of a symmetrical electrolyte (a") is defined as: 
a" = (a+ a. )' A = a 2 w 


The individual ionic activity coefficients are: y + =—±- 


m_, 


a 

y=— 

m. 


The mean ionic activity coefficient is: -y + =(y + -y _) 2 =( a +' a - \z — AL 




Ionic Strengh - is a useful generalisation for assessing the combined effects of the activities of several 
electrolytes in a solution of a given electrolyte. The value is given by the expression: 


/ = m i z ] Summation is over all ions in solution, both positive and negative. 

where nii is the molality of the i th ion and 
z is the charge on i th ion. 


6.7.3 Debye-Huckel Theory of Activity Coefficients 

Even in very dilute solutions of electrolytes there is long-range electrostatic force between charged ions, 
resulting in the lowering of the activity coefficients of the ions. These effects are evaluated by the Debye- 
Hiickel Theory. In dilute solutions Debye-Huckel expression states that the individual activity coefficient 
is given by: 

Az \fl where I - ionic strength 

— log -y.=-—-— z - charge on the ion 

1 + a, B V/ A, B - are constants characteristic of the solvent at a 

given temperature and pressure 
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a” - is a value which depends on the "effective diameter" 
of ions in solution and is determined by experiment. It is 
larger than values given for crystals because of solvation 
of ions 

When ionic strength becomes very small, — log y, = Az t V/ 

6.8 Solid Solutions and Melts 

-Review of atomic structure 
-Review of the structure of molecules 
-Types of bonds in crystals 
-Metallic crystals/crystals in general 
-Imperfections in crystals (Lattice Defects) 

(i) Point defects, ie. vacancies or holes from interstitial atoms, eg. IT in steel, 

(ii) Line defects, ie. dislocations 

(iii) Surface and interface defects, ie grain boundaries and crystallite boundaries. 

Solid solutions are common in metallic systems (alloys) and in inorganic systems, some at elevated 
temperatures. Some solid solutions have been reported in the literature about organic systems as well, 
even though these are not common. Depending on the locations of the foreign atoms (solutes) in the 
crystal lattice, solid solutions can either be substitutional (replacement of atoms of the parent metal or 
solvent) or interstitial (atoms occupy interstitial positions in the crystal) 

Two or more metals can form a homogeneous mixture in a single phase, ie. a solution. For example 10 
mass% Cu in Fe forms a homogeneous alloy. But as the Cu content exceeds 10% a Cu-rich phase is also 
formed in addition. In all these cases the lattice parameter (d) varies as more solute is added. 

Minerals are invariably solid solutions, eg. MnC0 3 (rhodochrosite) usually contains dissolved calcium, or 
Fe(II), Mg(II), or Zn(II) ions. The extent and kind of solution varies depending on the composition and 
crystal structure of the mineral acting as solvent. Most minerals are salts, crystalline compounds, made up 
of cations and anions at positions in the crystal structure. By the substitution of ions in the crystal 
structure solid solutions may be formed, eg. (Fe x Mni. x )C0 3 , or Fei_ x Al x OOH. 

The extent of solution and stability of the solution depends on the similarity or otherwise of the ions (eg. 
in ionic radius). With increase in temperature, as the thermal motion of the ions increase and with it the 
effective ionic radii, and therefore the effective ionic radii of all the elements become more nearly the 
same. The mutual solubility of minerals become greatly enhanced, eg. KC1 and NaCl are mutually soluble 
in all proportions above 500EC (r K + = 1.33A 0 , r Na+ = 0.95A°). 

As is known the activity of pure substances at standard state is 1, but that of (Fe x Mni_ x )C0 3 is not 1. For 
(Fe 0 . 2 Mno. 8 )C 0 3 , which might be considered as a dilute solution (y =1): 

a -x - °' 8 -Qp, 

u MnCO~ x MnCO~ Q 0+Q g — U '° 


Relatively few data are available on the behaviour of solid solutions and melts of geological interest. 
Flowever, those studied so far behave nearly as regular solutions, for which: 


logy,= 


B 


2.303 RT 


x\-B '(l-xj 2 


For most geological systems of interest, 0 < B'< ±5. In most cases the activity coefficient (y) approaches 
unity rapidly as the mole fraction approaches unity (1). Therefore even for solutions exhibiting large 
departures from ideal behaviour over much of the composition range the solvent may be considered to 



behave ideally as the solution becomes dilute, especially for minor substitutions in the mineral. It is 
reasonably safe to assume ideal behaviour if Xj$0.9. In the absence of the knowledge of B' no good 
predictions can be made if x, < 0.9. 

The activities of solids can be corrected as a result of compositional variations. For example: 

Problem : If calcite contains 5 mole% manganese carbonate, how would the product of the activities of 
calcium and carbonate ions change in aqueous solution in equilibrium with calcite at 25° C. 

Solution : 

CaC0 3s = Ca-Tq + C0 3 2 aq 

AG 0 = AG°fCa+ + AG°fC032- ■ AG°fCaC03 

= -132.18 + (-126.22) - (-269.78) = 11.38 
But AG° = -RT logK = -1.364 logK 

11.38 = -1.364 log ([Ca 2+ ] [C0 3 2 ]) / [CaC0 3 ] 

[Ca 2+ ] [C0 3 2 ] = 10' 8 - 34 ([CaC0 3 ] = 1 for pure calcite ) 

If Cao.95Mno.os)C0 3 , then [CaC0 3 ] — Xc a co 3 = 0.95 
Then [Ca 2+ ] [C0 3 2 ] = 10’ 8 34 x 0.95 = 10' 8 36 
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7 ELECTROCHEMISTRY 

7.1 Introduction 

There are generally two broad categories of chemical reactions; atom transfer reactios and electron 
transfer reactions. In atom transfer reactions, the chemical reaction does not lead to changes in the 
oxidation state of the constituent atoms of the reactants. However electron transfer reactions result in the 
change in the oxidation state of the reactants. These reactions may occur in homogeneous or 
heterogeneous systems. Electrochemistry is that branch of science that deals with chemical changes that 
are caused by electric current and chemical reactions that are accompanied by the transfer of charges 
(change of oxidation state of some or all the participating atoms). Such electron transfer reactions form 
the basis of some important industrial processes such as electroplating, electrowinning and 
electrorefinning, as well as some important geochemical processes that are involved in the formation of 
minerals or in determining the quality of groundwater. Other examples include the generation of the flow 
of electricity in car batteries as a result of chemical reactions in the cells, or the generation of electric 
current in solar cells, among others. It also helps to explain, for example the corrosion of metals. 

7.2 The Energetics of Electrochemical Reactions 

The thermodynamic driving force in any electrochemical reaction is related to the ease with which the 
change in oxidation state takes place, ie. the ease of electron transfer. The free energy decrease is a 
manifestation of the energetics of electron transfer. The driving force for the electron movement is the 
potential difference ()E) and the chemical reaction is driven by the thermodynamic force ()G). 

If both electrical and chemical driving forces are involved at the same time in electrochemical reactions, 
then there must be some kind of relationship between these two forces. We have already learnt that for a 
system undergoing a reversible process at constant temperature and pressure the change (decrease) in free 
energy of the system equals the maximum work (other than work of expansion) done by the system. That 
is: 

AG (T,p) = w max 

or in the form of a differential equation, -dG = *w max 

For example, the electrons moving in the external circuit of a galvanic cell such as a car battery, causes 
current to flow which is capable of doing work, such as starting a car. When an electric charge is 
transported across a potential difference work is done and the work done is equal to the product of the 
charge transported (Q, coulombs) and the electric potential difference (E). If the charge is transported 
reversibly, for example in a cell, then the work done is: 

Wmax “ Z F E 

where z is charge on ion, 

F is Faraday's constant, 96500 C) 

E is the electromotive force (EMF), which is the algebraic difference between the 
electrodes of a cell if no current is flowing through the cell. 

If you consider the chemical reaction (cell reaction): 

Metal = metal ion + electron , 

work is done by the system, releasing electrons to flow resulting in work being done. The free energy 
change for the reaction is: 



A G = p i0n - u metal = nL “ M-Laj ~ ^ In a i0? + FT In a mrtfl/ 


That is: 


E - 




Brne/a/ 


zF 


FT , a ion 

- In- 

zF 


metal 


This reduces to : 


E = — = E° 
zF 


RT ln a ™ 


1 metal 


In general for the reaction: aA + bB = cC + dD , (which is made up of two half-cell reactions, 

aA = cC + ze" a and bB + ze" c = dD) 

A G= cp c + du D — a\i A — b\x B from which 

This expression is known as the Nernst EEquation 

Taking a Daniel Cell as an example, the cell reaction: 

Zn + Cu 2+ = Zn 2+ + Cu 
is made of two half-cell reactions: 

Zn = Zn 2+ + 2e 
Cu 2+ + 2e = Cu 

At equilibrium, p Zll 2+(in solution) = p Z n 2 +(in Zn electrode) 

and p Zn = Mz„ 2 + + 2|Je Zn 

Pzn is the chemical potential of the atoms in the electrode, p /n2+ is the chemical potential of Zn ions in the 
ZnS0 4 solution and in the electrode, : e Zn is the chemical potential of the electrons in the Zn electrode. 

Similarly, 

Mcu 2 +(in solution) = p C u2+(in Cu electrode) and 

Mcu 2 + + 2p, Cu = pcu 




The equilibrium value of p,, for both reactions is determined by the activities of the sulphates in solution. 
In a reversible transfer of electrons, work done is: 

*w max = zF(E 2 - Ei) (or *w max = zF(E 2 - Ei) dn; dn is moles of electrons transferred ) 

The work done in the electron transfer at constant T and p involves decrease in free energy, therefore: 

*w max = +AG = +(p 2 - Pi) (where p 2 = p e Zn and p : = p e Cu ) 
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Combining the above equations gives: 

(jj 2 - Mi) = zF(E 2 - Ei) This expression gives the relationship between the chemical potential 

difference and the electrical potential difference for electron transfer. 

The relationship therefore between the chemical driving force and the electrical driving force therefore is: 

A G= — zFE 

If an external opposing electrical potential equal in magnitude to E and equal to the chemical driving 
force is applied across the electrodes of the cell then equilibrium is achieved. If this external voltage is 
decreased current flows in the external circuit from the Zn electrode to the Cu electrode and an equivalent 
ionic current (S0 4 2 ) flows through the cell at an equal rate. Zn 2+ ions are added by the corrosion of the Zn 
electrode to the electrolyte or Cu 2+ ions are removed by the deposition of Cu on the Cu electrode from the 
electrolyte. The reaction proceeds (CuS0 4 concentration decreases and ZnS0 4 concentration increases) 
causing the EMF to decrease until a new equilibrium is reached. 

In our earlier discussions we have indicated that some thermodynamic properties are extensive (such as 
U, S, N, etc), whilst others (such as temperature and pressure) are intensive. Potential difference, just as 
temperature, is an intrinsic property, which is independent of the quantities of the reacting species. We 
therefore simply cannot add potentials together when combining electrochemical reactions. Instead, the 
change in free energy, which is an extensive property and therefore depends on the quantities of the 
reactants can be calculated from the combination of the )Gs of the half-cell reaction. 

The effect of concentration on EMF 

The effect of concentration on the EMF can be seen through the Nernst equation, which is given in 
general for the cell reaction, 

aA + bB = cC + dD as E = E° - — 

zF 


ln- 


a r 


American and European Conventions 
With American convention: 

A G°= —zF<o 

With European convention: 

AG° = + zF%° = -zFE 

This means that E - - <o 


7.3 Electrochemical Cells and Electrode Potentials 
7.3.1 Daniel Cell 

The Daniel cell has already been discussed above. 



7.3.2 Galvanic cell 


It is a system capable of performing electrical work as a result of the occurence of chemical reaction, eg 
galvanic cell or car batteries (Pb plates between Pb0 2 sieves, all in sulphuric acid solution). 

At the anode: Pb(s) + HSCV = PbS0 4 (s) + H + + 2e 

At the cathode: PbOTsl + HSO/ + 3H + + 2e = PbSOTsI + 2H,Q 

Cell reaction: Pb(s) + Pb0 2 (s) + 2HS0 4 ' + 2H + = PbS0 4 (s) + 2H20 

The overall cell reaction indicates the consumption of acid and the production of water, leading to the 
change in density. 


7.3.3 Concetration cell 

It is a cell with identical electrodes inserted into solutions differing in only concentration. An example is 
the 02 concentration cell using lime-stabilised zirconia as a solid electrolyte: 

02(g, pressure I), Pt | Ca0-Zr02 | Pt, 02(g, pressure II) 

for which ^ _ _ RT ^ Pom) 

4F P 02[ll) 

The sudent is referred to discussions in the literature for other types of cells. 


7.3.4 Single electrode potentials 

When a peice of metal is immersed in an lectrolyte, the metal is called an electrode. The dissolution of the 
metal is represented by the cell reaction: 

Me = Me z+ + ze 

If the electrons are more “soluble” than the metal ion, then the charge on the electrode surface will be 
positive, that is the electrode/solution surface will acquire a positive galvanic potential. On the other hand 
if the metal ion is more soluble the electrode surface will be negatively charged. 

Electrodes are termed anodes or cathoides depending on the direction of flow of electrical charge or as to 
what whether the reaction taking place on the electrode is an oxidation (electron donation) or reduction 
(electron acceptance) reaction. When negative charge flows from the electrolyte to the electrode, an 
anodic current is said to be flowing and the electrode is said to be an anode (or when oxidation is taking 
place at the electrode). When electrons leave the electrode and pass into the electrolyte, a cathodic current 
is said to be flowing and the electrode is said to be a cathode (in other words, where reduction takes 
place). 


7.3.5 Standard Electrode Potentials 

It is impossible to construct a cell with a single electrode and measure its EMF. In order to measure the 
EMF, we choose a particular single standard electrode and assign to it a zero potential. Such a standard 
electrode is the standard H 2 electrode (H 2 electrode with gas at 1 bar pressure in contact with an aqueous 
solution of unit H + activity). The standard potential of the metal in comparison with the hydrogen 
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electrode, consists of pure metal in contact with an aqueous solution of the metal ion of unit activity. 
For example to measure the electrode potential of Zn the following cell is used. 

Zn | ZnCl 2 (a Z n 2 + = 1) || HCl(a H+ = 1) | H 2 (1 bar, Pt), 

for which the cell reaction is, 

Zn + 2H + = Zn 2+ + ViH 2 


E = E. 


anodic cell 


+ E, 


cathodic cell 


= E; 


anodic cell 


+0.762 V at 25EC 


(A list of standard potentials can be found in the Table of Electrochemical Series) 

7.4 Areas of the application of Electrochemistry. 

It is clear from the above that metals can dissolve in solutions or melts after giving out electrons if the 
electrochemical environment is favourable. The driving force is the EMF/potential difference and the free 
energy change()G). This knowledge is applied in numerous areas, including corrosion of metals and its 
prevention, the dissolution of metals or their compounds from minerals (leaching), or the production of 
metals or their compounds from solutions and melts (electrolysis, electrowinning and electrorefinning, 
electrodeposition), among others. 

7.5 Electron activity 

The electron activity is defined as: 

p, = -log (e) (compare with pH = -log {an + } ) 


The measure of electron activity gives an indication of whether an electrochemical reaction is taking 
place under an oxidising or a reducing condition. Usually when the value of the p, is large and positive 
then the conditions are strongly oxidizing. However, when the value of p, is small and negative then 
strong reducing conditions exist. 


For the redox reaction, 

aA + bB + ne = cC + dD, 


the activity quotient, Q is given by: 


[a c } c {aD\ d 

[a A ] a {a B } b \e\ 


At equilibrium Q e = K and log Q e = log K , therefore: 


log Q - log- 


a c \ -[a D \ 


[a A \ -[a B \ 


+ n pz 


pz = —log K — —log 


a c\ ■Va D \ 


p £ = p£° ~ — log - 


[a A \ -[a B \ 


Example : Calculate pe for a solution containing 0.00001 M Fe 3+ and 0.001 M Fe 2+ if t for the reaction: 
Fe 3+ + e = Fe 2+ the log K is 12.53. 



7.6 Electrolysis 


7.6.1 Faraday’s Laws 

The fundamentals of electrolysis have been laid down in FARADAY'S LAWS. 

1 st Law: The amount of chemical decomposition (ie. mass of liberated substance) is 

proportional to the quantity of electricity passing through the solution (or melt) 
of electrolyte. 

m ~ It (m ~ lit, if current varies) 

where m is the mass of substance liberated, g or kg 
I is the current, amps (A) 
t is time, s 

m = qlt, 

where q is the electrochemical equivalent, and is characteristic of the element and is defined as 
the mass in grams of an element liberated by a coulomb of electricity. 

In general, 


q = 


A 


F.z 


A is the atomic mass 

F is Faraday's constant (96520 coulombs) 
z is the charge transferred 


The electrochemical equivalent of some of the elements are: 


h 2 

1.044 x 10' 5 

Ag 

1.118 x 10' 3 

Cu 

3.29 x 10' 4 

Au 

6.805 x 10' 4 

A1 

9.319 x 10 5 

2 nd Law: 

If the same quantity 


electrolytes, the masses of substances liberated are proportional to their chemical 
equivalents: 


m 


m, 


We have seen earlier that when a metal is immersed in its salt solution and equilibrium chemical process 
takes place with the liberation of electrons: 

Me = Me z+ + ze 

The electrode potential for this reversible process can be thermodynamically derived for small currents. 
This equilibrium potential can also be measured using the calomel electrode, ie. 

~~ - (j> kal ( E e E m - Ekal ) 

When a large current passes across the boundary between the electrode and the solution, there is a shift in 
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the equilibrium potential and this shift resulting in a potential difference between these values is known 
as polarization. The difference between the equilibrium potential (h q ) and that of the polarised potential 
(E p ) gives the overpotential (overvoltage), r|: 

Ep - E e q — H 

Polarisation is due to the following: 

1) reaction(s) taking place at the electrodes (activation polarization). 

2) diffusion of ions from the body of electrolyte to areas around the electrode 
(concentration polarization). It can be reduced by agitation and movement of the 
electrolyte. 

3) insoluble films formed on the electrode causing ohmic resistance (resistance 
polarization). 

7.6.2 Minimum Electrolysing Voltage 

z z ^ j 

If you consider the reaction: M + ~+^H 2 0 = ~^0 2 + zH + + M 

which can be expressed as the sum of two half-cell reactions: 

M +: +ze~ = M and = -0 7 + ze~ L 

2 4 2 

the free energy change can be expressed as: A - ^[i Hi0 

The calculated free energy change for these reactions (see Table given), is positive and therefore does not 
favour spontaneous reaction in the direction of metal deposition. 


Metal ion 

AG° f (kJ/mol) 

Metal ion 

AG°t (kJ/mol) 

Ag + 

41.382 

Ni 2+ 

282.15 

Au 3+ 

84.436 

Pb 2+ 

260.832 

Cr 3+ 

151.316 

Pd 2+ 

105.754 

Cu 2+ 

171.38 

Sn 2+ 

263.758 

Fe 2+ 

328.966 

Zn 2+ 

383.724 


The ability to achieve metal deposition via electrolysis is based on the fact that the chemical potentials of 
the respective electrons associated with the cathodic and anodic reactions can be altered sufficiently to 
reverse the sign of the free energy change corresponding to the overall reaction. This can be done by 
applying an external source of electrical energy (eg from a DC supply such as a battery). 

By applying an external source for metal deposition to occur, the electrons at the cathode and anode 
would have acquired different chemical potentials and the overall reaction in this case can be expressed 
as: 

M +z +^H 2 0 + ze~ c L |0 2 + zH +1 + M + ze~ L 











and the free energy change for electrolysis will be: 

^G 2 = + z\i H . + u M - \x M ~ - Hi o + z\i t ' - z^ 

AG 2 =AG[ + z\x e — z\x e 

In the absence of externally applied voltage, the cathode and anode electrons have the same chemical 
potential and therefore the free energy change is equal. On the other in the presence of applied potential 
and at equilibrium, AG 2 = O and therefore, 

AGi -z ( \- \i e ) eq 

Re-writing this in the form of the relationship between the electrical potential and chemical potential 
gives: 

A G { -Z (0 e< - Q e ) eq 

AG —zFE where E min is the minimum electrolyzing potential. 

1 min 


7.7 pH AND Eh Measurements 

7.7.1 pH Measurement 

Definition: 


pH = -log a H + = -log [H+] (involves the activity of a single ion, a quantity that cannot be directly 

measured). 

pH is measured today using an EMF cell, that consists of a glass electrode dipping into the tests solution, 
together with a reference electrode (usually calomel, Hg/Hg 2 Cl 2 , electrode in saturated KC1 solution) to 
complete the circuit. 

The glass electrode consists of a bulb of special glass containing an acid solution, and an inner electrode 
of fixed voltage, usually Ag/AgCl, to conduct electrons reversibly into and out of the solution. When the 
bulb is immersed in a solution, an electric potential is developed between the inner and outer solutions 
that is proportional to the logarithm of [H + ] in the external solution. Thus the reference electrode has a 
constant potential, while that of the glass electrode varies with the hydrogen ion activity: 

Glass electrode, sol. x | KCl(sat), Hg 2 Cl 2 (s),Hg(l) 

The half-cell potential of a hydrogen ion sensitive glass electrode is: 


E'=C 


H 


+ 23 °l RT log[H +l ] 


n , ^ 2.303 RT , 

E x = c h - -p -log pH 


C H is a constant and differs from one electrode to the other, and depends on the specific composition of 
glass used, the standard solution in the bulb and the particular reference electrode in the circuit. 
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The overall potential of the cell neglecting the liquid junction potential is: 

2.303 RT 


F — F' — F — C + 

^ x ^ x ^ ref ^ H T 


F 


log[« l - E, 


ref 


(*) 


If solution sol.x is replaced by a standard solution, then: 


„ . 2.303 RT TT 

E s = C H + --- pH s - E r 


p IT—s —‘ref 

Subtracting (**) from (*) and rearranging: 

F[E X — E J This now serves as the basis for the modern definition of pH. 


pH x = pH s - 


2.303 RT 


Accuracy of pH measurements depends on the standards used, and also the error due to the liquid 
junction potential. If standard and test solutions are similar, then pH is more accurate. 

The Glass electrode 

General purpose glass electrodes are most efficient for pH 0-11. Above this, special "high pH" electrodes 
are necessary for accurate work. Glass electrodes are calibrated against at least 2, preferably 3 different 
buffers in order to check the slope of the line relating to the EMF of the glass electrode, which is: 

— 2.303RT 
F 

Therefore at constant T, 


dE, 

d pH 


(2.303 RT) 

p The value of this at 25°C is 0.05916 V/pH 


Calomel electrode 

This contains mercury in contact with mercurous chloride, which is in turn in contact with a standard 
solution of KC1. The half-cell reaction and potentials are: 


2Hg(l) + 2Cl(aq) = Hg 2 Cl 2 (c) + 2e 

Hg 2 Cl 2 \ 


J7I J7° ^ i 

Eh = E + ^jrln 


[H g f [cr 


Mercurous chloride and mercury have unit activity at standard state. Since KC1 solution is in equilibrium 
with solid KC1, chloride activity is constant and therefore, 

J7i r° j . RT , 1 

Eh = E + 2F' n F 


At 25°C, Eh is 0.244 V or more accurately 

Eh = 0.2444 - 0.00066(t-25) V (t = temp. °C) 










7 . 7.2 Eh Measurement 

Oxidation potential is measured with an electrode pair consisting of an inert electrode and a reference 
electrode (eg. calomel). The inert electrode used most is the bright platinum (or Au electrode) The role of 
the calomel electrode is to supply a known EMF and make electric connection with the system to be 
measured. The inert electrode acts as an electron acceptor or donor to the ions in the measured solution, 
that is if the potential of the half-cell containing the dissolved species is greater or less than that of the 
calomel reference electrode. 

Example : Calculate the oxidation potential of a Fe 2+ / Fe 3+ aqueous solution containing [Fe 2+ ] = 0.001 and 
[Fe 3+ ] = 0.01. 

The cell may be represented by: 

Pt(inert) I Fe aq 2+ (0.01), Fe aq 3+ (0.001) I Hg 2 Cl 2(c) , KCl aq , sat | Hg (l) 

Solution 

The half-cell reactions are: 

2Fe aq ‘ + — 2Fe aq 3+ + 2e EhF e2 +,Fe 3 + 

2Hg(l) + 2Cl-aq = Hg 2 Cl 2(c) + 2e EI^ e£ 

The overall cell reaction is: 

2Fe aq 2+ + Hg 2 Cl 2(c) = 2Fe aq 3+ + 2Hg (1) + 2Cb aq 
and FMFceii = Eh Fe2 +,Fe3+ - Eh re f (Eh ref = 0.244 V at 25°C) 


Eh F +2 +3 - Eh J-, +2 „ +3" t " 

Fe Fe Fe Fe 


0.0592 


log 


[Fe +3 ] 2 

[Fe* 1 ] 2 


Substituting the appropriate values into the above (Eh° F e 2 +,Fe 3 + = 0.771 V) gives Eh Fe2 +,Fe 3 + = 0.830 V. 
Therefore, 


EMF C eii = 0.830 - 0.244 V = 0.586 V 


(this would be measured if liquid junction 
potential is neglected). 


In general for an oxidation-reduction system: 
Re = Ox + ne 


Eh 


T7I ° I FF 1 

Eh + —log 
nr 


[Ox] 

[*] 


and 


EMFceii - EltRe.Ox - Eh re f 


( Eh Re ,ox is calculated from Nernst equation) 


Individual ionic activities can also be measured using Eh measurement and using their individual 
electrodes, eg. Ag-AgCl for Cl' ion, Ag-AgBr for Br", Ag-Ag 2 S for S 2 " ion activities, etc: 


O 


RT 

nF 


In a x ; 


(E° - Eh)-nF 

x = e x P 


Eh - E 


a 


RT 
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EMF measurement at high T and P : Glass electrodes can be used for EMF measurements up to 150 °C 
and 1000 bars pressure, whilst Ag-AgCl electrode can be used up to 250°C and 40 bars 



